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In a well-known investigation Sir William Thomson has discussed the problem of 
the bodily tides of a homogeneous elastic sphere, and has drawn therefrom very 
important conclusions as to the meat rigidity of the earth. ^"'' 

Now it appears improbable that the earth should be perfectly elastic ; for the con- 
tortions of geological strata show that the matter constituting the earth is somewhat 
plastic, at least near the surface. We know also that even the most refractory metals 
can be made to flow under the action of sufficiently great forces. 

Although Sir W. Thomson's investigation has gone far to overthrow the old idea of 
a semi-fluid interior to the earth, yet geologists are so strongly impressed by the fact 
that enormous masses of rock are being, and have been, poured out of volcanic vents in 
the earth's surface, that the belief is not yet extinct that we live on a thin shell over 
a sea of molten lava. Under these circumstances it appears to be of interest to inves- 
tigate the consequences which would arise from the supposition that the matter 
constituting the earth is of a viscous or imperfectly elastic nature ; for if the interior 
is constituted in this way, then the solid crust, unless very thick, cannot possess 
rigidity enough to repress the tidal surgings, and these hypotheses must give results 
fairly conformable to the reality. The hypothesis of imperfect elasticity will be prin- 

^ Sir William states tlia(i M. Lamib liad treated tlie subject at an earlier date, but in an entirely 
different manner. I am not aware, however, that M. Lamtii had fnlly discussed the subject in its physical 
aspect. 
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cipally interesting as showing how far Sir W. Thomson's results are modified by the 
supposition that the elasticity breaks down under continued stress. 

In this paper, then, I follow out these hypotheses, and it will be seen that the 
results are fully as hostile to the idea of any great mobility of the interior of the 
earth as is that of Sir W. Thomson". 

The only terrestrial evidence of the existence of a bodily tide in the earth would be 
that the ocean tides would be less in height than is indicated by theory. The subject 
of this paper is therefore intimately connected with the theory of the ocean tides. 

In the first part the equilibrium tide-theory is applied to estimate the reduction 
and alteration of phase of ocean tides as due to bodily tides, but that theory is 
acknowledged on all hands to be quite fallacious in its explanation of tides of short 
period. 

In the second part of this paper, therefore, I have considered the dynamical theory 
of tides in an equatorial canal running round a tidally-distorted nucleus, and the 
results are almost the same as those given by the equilibrium theory. 

The first two sections of the paper are occupied with the adaptation of Sir W. 
Thomsois3"'s work'"' to the present hypotheses ; as, of course, it was impossible to repro- 
duce the whole of his argument, I fear that the investigation will only be intelligible 
to those who are either already acquainted with that work^ or who are willing to 
accept my quotations therefrom as established. 

As some readers may like to know the results of this inquiry without going into 
the mathematics by which they are established, I have given in Part III. a summary 
of the whole, and have as far as possible relegated to that part of the paper the 
comments and conclusions to be drawn. I have tried, however, to give so much 
explanation in the body of the paper as will make it clear whither the argument is 
tending. 

The case of pure viscosity is considered first, because the analysis is somewhat 
simpler, and because the results will afterwards admit of an easy extension to the case 
of elastico-viscosity. 



I. 

THE BODILY TIDES OF VISCOUS AND ELASTICO-YISCOUS SPHEROIDS. 

1. Analogy hetiveen the flow of a viscous hody and the strain of an elastic one. 

The general equations of flow of a viscous fluid, when the effects of inertia are 
neglected, are 

* His paper will be found in Phil. Trans., 1863, p. 573, and §§ 733-737 and 834-846 of Thomson and 
Tab's 'Natural Philosophy,' edit, of 1867, 
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where x, y, % are the rectangular coordinates of a point of the fluid ; a, ^, y are the 
component velocities parallel to the axes ; p is the mean of the three pressures across 
planes perpendicular to the three axes respectively ; X, Y, Z are the component forces 

acting on the fluid, estimated per unit volunie ; v is the coefiicient of viscosity ; and 

d^ d d 

V^ is the Laplacian operation "t^+t^ + tt 
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Besides these we have the equation of continuity --—{—- — | — ~=0 

^ , '^ dx cly dz 

Also if P, Q, R, S, T, U are the normal and tangential stresses estimated in the 

usual way across three planes perpendicular of the axes 
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Now in an elastic solid, if a, /S, y be the displacements, m— J^^ be the coefficient of 
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dilatation, and n that of rigidity, and if S=y7+~p + 7^; the equations of equilibrium 
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P=:(m-7^)S+2fl^, Q=(m-n)S+2?^'^^, E=(m—n.)S+2«|j . . (4) 
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and S, T, U have the same forms as in (2), with n wiitten instead of v. 

"I '?2> 

Therefore if we put — p=5(P+Q+E')? we havej^9='— (m— -")S, so that (3) may be 
written 
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* Thomson and Tait's * Kat. PMl.,' § 698, eq. (7) and (8). 
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+7lV^a+X=0, &c.^ &c. 
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Also 



P:^ ^J_27^---, Q= &c., E;= &c. 

Tir dx 

Now if we suppose the elastic solid to be incompressible, so that m is infinitely 
large compared to n, then it is clear that the equations of equilibrium of the incom- 
pressible elastic solid assume exactly the same form as those of flow of the viscous 
fluid, n merely taking the place of v. 

Thus every problem in the equilibrium of an incompressible elastic solid has its 
counterpart in a problem touching the state of flow of an incompressible viscous fluid, 
when the effects of inertia are neglected ; and the solution of the one may be made 
applicable to the other by merely reading for ^'displacements" '^ velocities," and for 
the coefficient of '' rigidity " that of '' viscosity/^ 



2. A sphere under influence of bodily force. 



Sir W. Thomson has solved the following problem :- 



To find the displacement of every point of the substance of an elastic sphere exposed 
to no surface traction, but deformed infinitesimally by an equilibrating system of forces 
acting bodily through the interior. 

If for ''displacement" we read velocity, and for "elastic" viscous, we have the 
corresponding problem with respect to a viscous fluid, and mutatis mutandis the 
solution is the same. 

But we cannot find the tides of a viscous sphere by merely making the equilibrating 
system of forces equal to the tide-generating influence of the sun or moon, because the 
substance of the sphere must be supposed to have the power of gravitation. 

For suppose that at any time the equation to the free surface of the earth (as the 

viscous sphere may be called for brevity) is r=^a-{-^o-i, where cri is a surface harmonic. 

2 

Then the matter, positive or negative, filling the space represented by Scr^ exercises 
an attraction on every point of the interior ; and this attraction, together with that of 
a homogeneous sphere of radius a, must be added to the tide-generating influence to 
form the whole force in the interior of the sphere. Also it is a spheroid, and no 
longer a true sphere with which we have to deal. If, however, we cut a true sphere 
of radius a out of the spheroid (leaving out So"^), then by a proper choice of surface 
actions, the tidal problem may be reduced to finding the state of flow in a true sphere 
under the action of (i) an external tide-generating influence, (ii) the attraction of the 
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true sphere, and of the positive and negative matter fiUing the space Sc7i, but (iii) 
subject to certain surface forces. 

Since (i) and (ii) together constitute a bodily force, the problem only differs from 
that of Sir W. Thomson in the fact that there are forces acting on the surface of the 
sphere. 

Now as we are only going to consider small deviations from sphericity, these surface 
actions will be of small amount, and an approximation will be permissible. 

It is clear that rigorously there is tangential action ^'^ between the layer of matter 
XcTi and the true sphere, but by far the larger part of the action is normal, and is 
simply the weight (either positive or negative) of the matter which lies above or 
below any point on the surface of the true sphere. 

Thus, in order to reduce the earth to sphericity, the appropriate surface action is 
a normal traction equal to -^gwXcTi, where g is gravity at the surface, and lu is the 
mass per unit volume of the matter constituting the earth. 

In order to show what alteration this normal surface traction will make in 
Sir W. Thomson's solution, I must now give a short account of his method of 
attacking the problem. 

He first shows that, where there is a potential function, the solution of the problem 
may be subdivided, and that the complete values of a, /3, y consist of the sums of two 
parts which are to be found in different ways. The first part consists of any values of 
a, ^, y, which satisfy the equations throughout the sphere, without reference to surface 
conditions. As far as regards the second part, the bodily force is deemed to be non- 
existent and is replaced by certain surface actions, so calculated as to counteract the 
surface actions which correspond to the values of a, /3, y found in the first part of the 
solution. Thus the first part satisfies the condition that there is a bodily force, and 
the second adds the condition that the surface forces are zero. The first part of the 
solution is easily found, and for the second part Sir W. Thomson discusses the case of 
an elastic sphere under the action of any surface tractions, but without any bodily 
force acting on it. The component surface tractions parallel to the three axes, in this 
problem, are supposed to be expanded in a series of surface harmonics ; and the 
harmonic terms of any order are shown to have an effect on the displacements inde- 
pendent of those of every other order. Thus it is only necessary to consider th^ 
typical component surface tractions Ai, B^, C^- of the order i. 

He proves 'that (for an incompressible elastic solid fqr which m is infinite) this one 
surface traction A^, B^, C^ produces a displacement throughout the sphere given by 
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(5)t 



* I shall consider some of the effects of this tangential action in a future paper, viz. : " Problems con- 
nected with the Tides of a Viscous Spheroid," read before the Rojal Society qn December 19th, 1878. 
t Thomson and Tait's 'Nat. Phil.,' 1867, § 737, equation (58). 
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with symmetrical expressions for /3 and y ; where '^ a,nd ^ are auxiliary functions 
defined by 



^/ 
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In the case considered by Sir W. Thomson of an elastic sphere deformed by bodily 
stress and subject to no surface actiori, we haye to substitute in (5) and (6) only those 
surface actions which are equal and opposite to the surface forces corresponding to the 
first part of the solution f' but in the case which we now wish to consider, we must 
add to these latter thQ components of the normal traction —gwt(f%, and besides must 
include in the bodily force both the external disturbing foi-ce, and the attraction of 
the matter of the spheroid on itself. 

Now from, the forms of (5) and (6) it is obvious that the tractions which correspond 
to the first part of the solution, and the traction -~g^t{;So-i produce quite independent 
eflPects, and therefore we need only add to the complete solution of Sir W. Thomson's 
problem of the elastic sphere, the terms which arise from the normal traction -r-gwXcri. 
Finally we must pass from the elastic problem to the viscous one, by reading v for n, 
and velocities for displacements. 

I proceed then to find the state of internal flow in the viscous sphere, which results 
from a normal traction at every point of the surface of the sphere, given by the 
surface harmonic S^-. 

In order to use the formulae (5) and (6), it is first necessary to express the 

component tractions - S^-, '- S,-, - S,- as siux^face harmonics, 

-'- Cl> (If (m 

Now if Yi be a solid harmonic, 
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dx 
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The quantities within the brackets [ ] being independent of r^ and being 
surface harmonics of orders i—l and ^+1 respectively, we have - S^- expressed as the 

Co 

sum of two surface harmonics A^-^i, -^i+i, where 



* Where tlie solid is incompressible, this surface tractioipL is normal to the sphere at every point, 
provided that the potential of the bodily force is expressible in a series, of solid harmonics. 
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^2-1 O' 
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Similarly - Si, - Sj may be expressed as B2«.]^+Bi+i Q'^id Q^i+Q+i? where the B's 

Cv Clr 

and C s only differ from, the As in having y, z written for x. 

We have now to form the auxiliary f mictions "^i^i, #i corresponding to A^-^^j B^^j, 
Q^i and ^i, $,-+3 corresponding to Af^^, B^+i, Q+x- 

Then by the formulae (6) 
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Thus 
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Then by (5) we form a corresponding to A^-^^j^, B^^^, Q^x, and also to Ai^^, ^i-hif 
Q^25 and add them together. The final result is that a normal traction Si gives, 
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(7) 



and symmetrical expressions for yS^ and y^ 

a'5 ^', y are here written for % /J^ y to show that this is only a partial solution, 
and V is written for n to show that it corresponds to the viscous problem. If we 
now put Si= —gwcTi, we get the state of flow of the fluid due to the transmitted 
pressure of the deficiencies and excesses of matter below and above the true spherical 
surface. This constitutes the solution as far as it depends on (iii). 

There remain the parts dependent on (i) and (ii), which may for the present be 
classified together ; and for this part Sir W. Thomson's solution is directly applicable. 
The state of internal strain of an elastic sphere, subject to no surface action, but 
under the influence of a bodily force of which the potential is W^-, may be at once 
adapted to give the state of flow of a viscous sphere under like conditions. The 
solution is — 
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with symmetrical expressions for /3'' and y\ 

I will first consider (ii) ; i.e., the naatter of the earth is now supposed to possess the 
power of gravitation. 

The gravitation potential of the spheroid r=a+cr^ (taking only a typical term 
of cr) at a point in the interior^ estimated per unit volume, is 






2a 



2^^-l\a 



cr; 



according to the usual formula in the theory of the potential. 

Now the first term, being symmetrical round the centre of the sphere, can clearly 
cause no flow in the incompressible viscous sphere. We are therefore left with 



cr,:. 



Now if — — -( -j cr/ be substituted for W,- in (8), and if the resulting expression be 

3 
compared with (7) when -—giocri is written for S,-, it will be seen that — a''= r-^ — ~a\ 



Thus 
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• (9) 



with symmetrical expressions for yS'+/3'' and y-^-y'. 

Equation (9) then embodies the solution as far as it depends on (ii) and (iii). And 

2 . . . . 

since (9) is the same as (8) when— -('?*— 1)V« is written for W,-, we may include all 

the effects of mutual gravitation in producing a state of flow in the viscous sphere, by 
adopting Thomson's solution (8), and taking instead of the true potential of the layer 

* 'Nat. Pliil.', § 834, equation (8) wlien m is infinite compared -with n, and ^ — 1 written for ?', and v 
replaces n. 

t The case of § 815 in Thomson and Tait's 'Nat. Phil.' is a special case of this. 
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2 . ■ . 
of matter cr^, — -({— l) times that potential, and by adding to it the external disturbing 

potential. 

We have now learnt how to include the surface action in the potential ; and if W^ 
be the potential of the external disturbing influence, the effective potential per unit 
volume at a point within the sphere, now free of surface action and of mutual 

gravitation, is W^ — ^ ~ (r^-=r'Ti suppose. 

The complete solution of our problem is then found by writing r^T^ in place of W^ 
in Thomson's solution (8).'''" 

In order however to apply the solution to the case of the earth, it will be convenient 
to use polar coordinates. For this purpose, write wr^S^ for W^, and let r be the radius 
vector ; 6 the colatitude ; ^ the longitude. Let p, w, v be the velocities radially, and 
along and perpendicular to the meridian respectively. Then the expressions for p, ot, v 

d d 

will be precisely the same as those for a, ^, y in (8), save that for — we must put — ; 

f^ Cv vl' -^ p Cv Cv 

ay rsm ad^ dz rdd 

Then after some reductions we have 



P~ 2{i-l)[2(i+iy + l]v ' 

^~ 2{i-l)[2(i + iy + l]v ''"' dd 

i{i + 2)a^-{i-l)(i + 2,y^ r'-i dTi 
' 2{i-l){2(i+lf+l'\v sin dd4> 



S 



V. 



(10)t 



where Ti=wi S,-— 20777-7-7 " • 

These equations for p, t*t, v give us the state of internal flow corresponding to the 
external disturbing potential r'S,-, including the effects of the mutual gravitation of 
the matter constituting the spheroid. 

* The introduction of tlie effects of gravitation may be also carried out syntlietically, as is done by Sir 
W, Thomson (§ 840, 'Nat. Phil.') ; but the effects of the lagging of the tide-wave render this method 
somewhat artificial, and I prefer to exhibit the proof in the manner here given. Conversely, the elastic 
problem may be solved as in the text. 

t There seems to be a misprint as to the signs of the (^'s in the second and third of equations (13) of 
§ 834 of the 'Fat. Phil.' (1867). When this is corrected ju, and j/ admit of reduction to tolerably simple 
forms. It appears to me also that the differentiation of p in (15) is incorrect; and this falsifies the 
argument in three following lines. The correction is not, however, in any way important. 



JM.130Cyv/XjA-AX-A,. 
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3. The form of the free surface at any time. 



If p' be the surface value of p, then 



P = 



i(2i + l) 



w 



i+i 



2{i-l){2(i-rlf-¥l] V 



T. 



Hence after a short interval of time 8^, the equation to the bounding surface of 
the spheroid becomes r=:=a+o"^+p'S?^ ; but during this same interval^ o-^ has become 

0^5 whence 



dt 



or 
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% 



gwa 



2 (i+ 1)^ + 1 V 



(Ti 



% 



gioa 



-cn = 
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dt ' 2(^ + 1)^ + 1 V ' 2(^-l)[2(^+l)Hl] V 



kj^ • • • • \x.L/ 



This differential equation gives the manner in which the surface changes/ under the 
influence of the external potential r'S,-. 

If Si be not a function of the time, and if ^^ be the value of cxi when t=^0, 



0-; = 



2^ + 1 a% 
2(i-l) T 



1—exp 



-gwait 
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(12) 



* 



When t is infinite 



(Ti = 



2i + l a% 



2(^-1) g 



• • • * I 1 o I 



and there is no further state of flow, for the fluid has assumed the form wh^ich it 
would have done if it had not been viscous. This result is of course in accordance 
with the equilibrium theory of tides. 

If Si be zero, the equation shows how the inequalities on the surface of a viscous 
globe would gradually subside under the influence of simple gravity. We see how 
much more slowly the change takes place if i be large ; that is to say, inequalities of 
small extent die out much more slowly than wide-spread inequalities. Is it not 
possible that this solution may throw some light on the laws of geological subsidence 
and upheaval 1 

4. Digression on the adjustments of the earth to a form of equilibrium'. 

In a former paper I had occasion to refer to some points touching the precession of 
a viscous spheroid, and to consider its rate of adjustment to a new form of equilibrium. 



^ I write " exp." for " e to tlie power of." 
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when its axis of rotation had come to depart from its axis of symmetry. ^^ I propose 
then to discuss the subject shortly, and to establish the law which was there assumed. 
Suppose that the earth is rotating with an angular velocity o) about the axis of z, 
but that at the instant at which we commence our consideration the axis of symmetry 
is inclined to the axis of z at an angle a in the plane of a;^. and that at that instant 
the equation to the free su.rface is 



, ^ , 5m /I 



- — [cos a cos ^+sin a sin 6 cos <^]^ j V 



4 \ 

where m is the ratio of centrifugal force at the equator to pure gravity, and therefore 
equal to 



C()% 



9 



5ma /I 



Then putting {=2 in (12), and dropping the suffixes of S, s, cr, s-—~ i 
We may conceive the earth to be at rest, if we apply a potential 



l-^ ?)• 



so 



that 



wr^S=~o)hvrH ~— cos^ d 



S=-cu^f-— cos^ 6 



By (12) we have 



5a^S 



cr 



^9 
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Then, substituting for S and 5, and putting k= 



+s exp ( — 



2wga 
~1^ 



2wgat 



(7= -J- 4 (-— cos^ ^|[1— 6.3y(— Aci^)] + (o"~[<^^^^ a COS (9+sin a sin Q cos ^\^X'p{—Kt) 

Now 

[1 — ^^^jp(-~^0] ^^^ 6-\-exp{—Kt){coB a cos ^+sin a sin 6 cos ^Y 
= cos^ 6 [1 — sin^ a 6ccp(— fc^)]+sin^ a sin^ 9 cos^ ^ exply—Kt) 
+ 2 sin a cos a sin 6 cos 9 cos j> exp{—'Kt). 

Therefore the Cartesian equation to the spheroid at the time t is, 



x^ + 'if + ^^ 



1 Hr 



5m 

T 



a 



2 



5 m 



— — {2;^(l — sin^ a exp{—Kt))-\-x^ sin^ a exp{—Kt)-^2xz sin a cos aexp[—Kt)] 



or 



^ " On the Influence of Geological Changes on the Earth^s Axis of Rotation," PhiL Trans., Vol. 167, 
rart 1., sec. o. 
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\'\'\.w^ aex'pi^—Kt) ^+j!/^+^^i 1 +— (1— sin^ae.T;p(— /c^)) \ 
+ 5m sin a cos a xz exp{ — Kt) = a^( 1 + "~ 

Let a' be the inclination of the principal axis at this time to the axis of z, then 

, sin 2a exp(—Kt) 

tan ZCL — 7~~~Tt '- 9 7 In 

1—2 snr a exfi^—tct) 

If a be small, as it was in the case I considered in my former paper, then 

a'=a expi — kI) and -7—== — Ka , 
■^ ^ at 

Therefore the velocity of approach of the principal axis to the axis of rotation varies 
as the angle between them, which is the law assumed. 

Also K=-—~~, so that K (the p of my former paper) varies inversely as the coefficient 

of viscosity, — as was a^bo assumed. 



5. Bodily tides in a viscous earth. 






The only case of interest in which S^- of equation (11) is a function of the time, is 
where it is a surface harmonic of the second order, and is periodic in time ; for this 
will give the solution of the tidal problem. Since, moreover, we are only interested 
in the case where the motion has attained a permanently periodic character, the 
exponential terms in the solution of (11) may be set aside. 

Let 82=8 cos {vt'{-7]), and in accordance with Thomson's notation,t let — =g;, and 
r — ^=r; and therefore ^— — =- 

Then putting i^= 3 in (11), and omitting the suffice of cr for bi'evity, we have 

-^+~^<^=-^^os{vt-j-7]) . (14) 

It is evident that cr must be of the form A cos (t't+B), and therefore 

A{— ^r sin ('z;?^+B) + gcos ('z;i^+B)}=aS cos {vt-\~7)) 

* In certain cases the forces do not form a rigoronslj equilibrating sjstem, but there is a very small 
couple tending to turn the earth. The effects of this unbalanced couple, which varies as the square 

> 771 

of o '^7 will be considered in a succeeding paper on the " Precession of a Viscous Spheroid." (Read 

before the Royal Society, December 19th, 1878.) 
t^Kat. Phil.', § 840, eq. (27). 
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or if we put taB e == — , 



Ag sec e cos ('y^+B+e) = aS cos {vt-r{-r]). 



ac 



Hence A=-S cos e, and B=')? — e, 
Therefore the solution of (14) is^ 



a, 



0"=:-S C0S6C0S ('y^-f-T; — e) ,,...,,. (15) 

<w— w<w><M £/ ♦f -4- ty %y t/ 

Where tan e=—=~ 

g Agaw 

But if the globe were a perfect fluid, and if the equilibrium theory of tides were 
true, we should have by (13), 

o"=--.aS COS {vt'{-7))=---G0B {vt-\-rj). 

Thus we see that the tides of the viscous sphere are to the equilibrium tides of a 

(• 
fluid sphere as cos e : 1, and that there is a retardation in time of -. 

A parallel investigation will be applicable to the general case where the disturbing 
potential is wr^ S^* cos {vt-\'rj) ; and the same solution will be found to hold save that 

we now have tan €= ^ , and that in place of fl we have -j~. — tt^. 

^ gaw ^ ^ (2^-M)<x 



6. Diminution of ocean tides on equilihrium theory. 

Suppose now that there is a shallow ocean on the viscous nucleus, and let us find 
the effects on the ocean tides of the motion of the nucleus according to the equilibrium 
theory, neglecting the gravitation of the water. 

The potential at a point outside the nucleus is 

S'7 + 59'(7J ^+^^^S cos (^^+7?), 

and if this be put equal to a constant, we get the form which the ocean must assume. 
Let r=^a+u be the equation to the surface of the ocean. Then substituting for r in 
the potential, and neglecting u in the small terms, and equating the whole to a 

constant, we find, 

3 
— c7tfc+-(7cr+a^S cos {vt'\-rj)=-Q 



or 



3 . a^ 



't^zz:-cr+— S cos {vt-^-ij), 
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But the rise and fall of the tide relative to the nucleus is given by it-— a, and 

u — (jz=L — COS (t'^+7?)—- -cr 

g ^ '^ b 

==- — [cos [vt-\-r})— cos €COS (1;^+'^"-"^)] 

:^— - — sin € sin ('^;?^+')7 — e) ..,....• (17) 

E'ow if the nucleus had been rigid, the rise and fall would have been given by 

" — cos (t'^+'^) = B[ cos (vt-\-rj) suppose. 

Therefore 

t/, — cr== — Hsin € sin (t'^+'>? — e) • (18) 

Hence the apparent tides on the yielding nucleus are equal to the tides on a 
rigid nucleus reduced in the proportion sin e : 1 ; and since •— sin ('^2^ +->;—€) 

= cos ('?;^+77+-— e) they are retarded by -(e— ~j. As e is necessarily less than -, 
this is equivalent to an acceleration of the time of high water equal to -( - — ej. 

It is^ however, worthy of notice that this is only an acceleration of phase relatively 

> sin 6^ po^ € 

to the nucleus, and there is an absolute retardation of phase equal to arc -tan ^—-^^ ^, 

O -T" Z COS 6 

7. Semidiurnal and fortnightly tides. 

Let the axis of z be the earth's axis of rotation, and let the plane of xz be fixed 

< 

in the earth ; let c be the moon's distance, and m its mass. 

Suppose the moon to move in the equator with an angular velocity o) relatively to 
the earth, and^let the moon's terrestrial longitude, measured from the plane of xz^ at 
the time t be o)t. 

Then at the time t, the gravitation potential of the tide generating force, estimated 
per unit volume of the earth's mass is 

— ■- ~-^W7'^\ -— sin^ 6 cos^ {(l)—o}t) 
which is equal to 



6 771 ^ / 1 ^ ^\ . 3 m- 



o 



, .wrH-— cos^ 0]+~ -wr^{w^^ 9 cos 26 cos 2o)t-\- sin^ 6 sin '26 sin 2(Dt}. 

The first term of this expression is independent of the time, and therefore produ.ces 
an effect on the viscous earth, which will have died out when the motion has become 
steady ; its only effect is slightly to increase the ellipticity of the earth's surface. 



diurnal tide is 7^, where tan e=-^ — = -----. Also the height of the tide is less than the 
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The two latter terms give rise to two tides, in one of which (according to previous 
notation) 

) oil 

S cos (vt-\-ri)^='-. -^ sin^ cos 2(^ cos 2a)t^ 
and in the second of which 

S COS {vt'\'7])=-'—~ -^sin^ 6^ sin 2</)C0s( 2a)i^+'^)- 

Now e, which depends on the frequency of the tide generating potential, will clearly 
be the same for both these tides ; and therefore they will each be equal to the corres- 
ponding tides of a fluid spheroid, reduced by the same amount and subject to the 
same retardation. They may therefore be recompounded into a single tide ; and 
since v will here be equal to 2a), it follows that the retardation of the bodily semi- 

e T , 2oyX 19u&) 

I—, where tan e=-^ — = 

2a> g gaio 

corresponding equilibrium tide of a fluid spheroid in the proportion of cos e to unity. 

Similarly by section (6) the height of the ocean tide on the yielding nucleus is given 

by the corresponding tide on a rigid nucleus multiplied by sin e, and there is an accelera- 

TT 6 

tion of relative high water equal to ~ -— . 

^ ^ 41:0) 2a) 

The case of the fortnightly tide is somewhat simpler. 

If O be the moon^s orbital angular velocity, and I the inclination of the plane of the 
orbit to the earth's equator, then the part of the tide generating potential, on which 
the fortnightly tide depends, is — 

'- ^ wr^ sin^ I ( 3 "~ cos^ 6) cos 2D.t 

and we see at once by sections (5) and (6) that tan e = — — . The bodily tide is the 

•^ ^ ' ^ ^ gaw *^ 

tide of a fluid spheroid multiplied by cos e ; the reduction of ocean tide is given by sin e ; 

and there is a time-acceleration of relative hisfh water of -777 — —77 or i- ^ of a 

week. 

In order to make the meaning of the previous analytical results clearer, I have 

formed the following numerical tables, to show the eflPects of this hypothesis on the 

semidiurnal and fortnightly tides. The coefiicient of viscosity is usually expressed in 

1 y 7jct) 

cavitation units of force so that the formula for e becomes, tan e = . In the 

^ • vja 

tables V is expressed in the centimetre-gramme-second system, and in gravitation uinits 
of force ; a is taken as 6*37 X 10^, and w as 5*5, and the angular velocity co of the moon 
relatively to the earth as '00007025 radians per second. 

With these data I find v= 10^^ X 2*625 tan e. As a standard of comparison with 
the coefficients of viscosity given in the tables, I may mention that, according to some 
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rough experiments of my own, the viscosity of British pitch at near the freezing 
temperature (34° Fahr.), when it is hard and brittle, is about 10^ X 1*3 when measured 
in the same luiits. 





Lunar 


Semidiurnal Tide. 








Height of 


Height of 


High tide 


Coefficient 


Betardation 


bodily tide is 


ocean tide is 


relatively to 


of viscosity 


of bodily tide 


1/ 

tide of 


tide on 


viscous nucleus 


xlO-^^ 




fluid spheroid 


rigid nucleus 


accelerated by 


(v x 10-^^) 


multiplied by 


multiplied by 




\ * 




(cos t). 


(sin £}. 




Hrs. min. 






Hrs. min. 


Fluid 





1^000 


•000 


3 6 


46 


21 


•985 


•174 


2 46 


96 


41 


•940 


•342 


2 25 


162 


1 2 


•866 


•500 


2 4 


220 


1 23 


'7 m 


•643 


1 44 


313 


1 44 


•643 


'766 


1 23 


465- 


2 4 


•500 


•866 


1 2 


721 


2 25 


•342 


•940 


41 


1,488 


2 46 


•174 


•985 


21 


Higid 00 


3 6 


•000 


1-000 





Jj'ortnightly Tide. 




Days. hrs. 






Days. hrs. 


Elnid 





I'OOO 


•000 


3 10 


1,200 


9 


•985 


•174 


3 1 


2,500 


18 


•940 


•342 


2 16 


4,000 


1 3 


•866 


•500 


2 6 


5,800 


1 12 


'766 


•643 


1 21 


8,300 


1 21 


•643 


'766 


1 12 


12,000 


2 6 


•500 


•866 


1 3 


19,000 


2 16 


•342 


•940 


18 


39,300 


3 1 


•174 


•985 


9 


Rigid 00 


3 10 


•000 


1^000 






I now pass on to a case which is intermediate between the hypothesis of Sir W. 
Thomson and that j ust treated. 

8. The tides of an elastico-viscous spheroid. 

The term elastico-viscous is used to denote that the stresses requisite to maintain 
the body in a given strained configuration decrease the longer the body is thus con- 
strained, and this is undoubtedly the case with many solids. In the particular case 
which is here treated, it is assumed that the stresses diminish in geometrical progres- 
sion, as the time increases in arithmetical progression. If, for example, a cubical 
block of the substance be strained to a given amount by a shearing stress T, and 

maintained in that position, then after a time t, the shearing stress, is Tej^pf— - 
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The time t measures the rate at which, the stress falls off, and is called (I believe by 
Professor Maxwell) " the modulus of the time of relaxation of rigidity ; ^' it is the 
time in which the initial stress has been reduced to e~^ or '3679 of its initial value. I 
do not suppose, however, that any solid conforms exactly to this law ; but I conceive 
that it is often useful in physical problems to discuss mathematically an ideal case, 
which presents a sufficiently marked likeness to the reality, where we are unable to 
determine exactly what that reality may be. 

Mr. J. G. Butcher has found the equations of motion of such an ideal sub- 
stance from the consideration that the elasticity of groups of molecules is continually 
breaking down, and that the groups rearrange themselves afterwards. ""''^ These con- 
siderations lead him to the following results for the stresses across rectangular planes 
at any point in the interior, viz. (with the notation of § 1) : — 



P=(„-„)a+.„('i+|)-'(^-+ i»), s=n(i+i)-(f +|; 



and similar expressions for Q, R, T, U ; where m— - is the coefficient of dilatation, n 

that of rigidity, S the dilatation, and a, ^, y, the components of flow. 

These expressions are clearly in accordance with the above definition of elastico- 

viscosity, for -+-=,^^~+-j. 

If the expressions for P, S, &c., be substituted in the equations of equilibrium 
of the elementary parallelopiped, it is found by aid of the equation of continuity 

dS da , dS , dy ^. ' . .... 1,1 

-—= - — L—- j--^^ that when mertia is nee^iected 

dt dx^ dy^ di ^ 



1 ^\-i 
t^ "^ dt. 



i^^-^*^^)l+^|JI+^V'4 "^ ^"^ 



and two similar equations. 

By the same reasoning as in § 1, we may put, ^^=-~~~;;^'> ^iid. the equations become 



1 ^\~i 
i^ dt) 



1 m d 

I 

t m — 



m—^ndt 



^-nV^aT+X=0. 

ax J 



Then supposing the substance to be incompressible, so that m is infinitely large com- 
pared to n^ and therefore m-^m—^n is unity, the equations become 



^^+,J1+|3 V^a+X=0 



dx 



dt. 



and two similar equations. 



* Proc. Loud. Math, Soc, Dec. 14, 1876, p. 107-9. It seems to me that the hypothesis ought to repre- 
sent the elastico- viscosity of ice very closely. 

MBCCCLXXIX. D 
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Now these equations have exactly the same form as those for the motion of a viscous 
fluid, save that the coefficient of viscosity v is replaced hjn (t+tJ • We may there- 
fore at once pass to the differential equation (11) which gives the form of the surface 
of the spheroid at any time. 

Substituting, therefore, in (11) for -, - (7+ i:)) we get 



1 + 



% gwa 



d(Ti i gwa i{2i-\-l) locd'^'^ /I d\^ 

In/' . iNO , 1 ~~7 ^i r>/.' i\rfi/.' . -IN 9. , i 1 ~~r \ ~Z\~7j. i^i* 



2(^ + l)Hl n \dt'^2{i-^iy-\-l ni^' 2(^-l)[2(^ + l)Hl] ^ \t"^ 

This eqiiation admits of solution just in the same way that equation (11) was solved ; 
but I shall confine myself to the case of the tidal problem, where i=-2 and 
82=8 cos (^'^+7;). In this special case the equation becomes 



2gwa\dcr 2gwa 5wa^ 



19^2./^^ ' 19nt 19n 



- cos {vt-\-rjD—v^m {vt-\~7}) 



S. 



And if we put h 1 =7, "ta.n xlj=vt, and 5=— 

^ zgwa k oa 

This may be written 

da , h vah c\ / , . i\ 

— +-0-= — : — r!S COS [vt-\-7)-±\b). 

In the solution appropriate to the tidal problem, we may omit the exponential term, 
and assume a=A cos {vt+B). Then if we put tan x= r 

—+-(7=-: — cos (t;?!+B+x). 
dt t sm Y ^ ^' 



X 



Whe^ it follows that 6=07+'/'— "X' ^^^ 



so that 



. a -, sm y a cos y 
g sm i/r 3 cos y 



a^ cos y , , , , \ 
cr= — ^ cos('?;i^+77+T//— v). 



Hence the bodily tide of the elastico- viscous spheroid is equal to the equihbrium tide 

COS 'V 

of a fluid spheroid multiplied by — ^, and high tide is retarded by x— ^-r-'?^. 

The formula for tan x niay be expressed in a somewhat more convenient form ; we 

have tan \b=v\ and therefore tan x=tan t//+— — 

Zigivcc 
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But nt is the coefficient of viscosity, and in treating the tides of the purely viscous 
spheroid we put tan €=r X coefficient of viscosity; therefore adopting the same 

notation here, we have tan ^=:tan i//+tan e. 

If the modulus of relaxation t be zero, whilst the coefficient of rigidity n becomes 
infinite, but nt finite, the substance is purely viscous, and we have i//=0 and x^=-^y ^^ 
that the solution reduces to the case already considered. If t be infinite, the sub- 
stance is purely elastic, and we have ti'=:7, y=o ^^.d since — ^znh- — Vj therefore 
^ ^ ^2^2 cos y^r sm i/r' 

a-=— S cos {vt'\'ri). 

But according to Thomson's notation'"^ =-, so that <t= S cos ('i^^+'>?X which is 

the solution of Thomson's problem of the purely elastic spheroid. 

The present solution embraces, therefore, both the case considered by him, and that 
of the viscous spheroid. 



9. Ocean tides on an elastico-viscous nucleus. 

If r=a-{-u be the equation to the ocean spheroid, we have, as in sec. (6), that the 
height of tide relatively to the nucleus is given by 

a^ 2 

U--or=—SGOs{vt-^r])---or, 

and substituting the present value of cr, 



Zi (Aj ^ 

u — cr=-- b 



cos X 



cos (^^+7?) - -rf ^^^ (^^+^+^-x) 



5 s L ^ ' cos i/r 

If the nucleus had been rigid the rise and fall would have been given by H cos (^'f +??), 

2 (X 
where H=- ~ S ; therefore on the yielding nucleus it is given by 

5 3 

COS y ^ ' ^^' 

■=. — H cos X (tan;5(— tan \//) sin (^'^+'>?— x) 
= — H cos X tan € sin ('I^^+t;-— x)- 

* ' Nat. Phil.', § 840. 
D 2 
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Hence the apparent tides on the yielding nucleus are equal to the corresponding tides 
on a rigid nucleus reduced in the proportion of cos x "tan e to unity, and there is an 

acceleration of the time of high water equal to -(- — x)' 

As these analytical results present no clear meaning to the mind, I have compiled 
the following tables. I take the two cases considered by Sir W. Thomson, where 
the spheroid has the rigidity of glass, and that of iron, and I work out the results for 
various times of relaxation of rigidity, for the semidiurnal and fortnightly tides. The 
last line in each division of each table is Thomson's result. 

I may remind the reader that the modulus of relaxation of rigidity is the time 
in which the stress requisite to retain the body in its strained configuration falls 
to '368 of its initial value. 

Spheroid with Rigidity of Glass (2-44: X 10^). 



TiUnar Semidiurnal Tide. 










High tide 


Modulus of 
relaxation of 


Coefficient of 

• ^ ' A 


Ocean tide 
is tide on rigid 


relatiyely to 
nucleus is 


rigidity 
(t). 


Viscosity 


nucleus 
multiplied by 


accelerated by 
\2 ~ ^j^ 






(cos X tan e). 


Hrs. 






Hrs. min. 


Fluid 





•000 


3 6 


1 


88 


•256 


1 44 


2 


176 


•342 


1 3 


3 


264 


•370 


45 


4 


351 


•382 


34 


5 


439 


•388 


28 


Elastic 00 


00 


•398 





Fortnightly Tide. 


Days, hrs. 






Days. hrs. 


Fluid 





•000 


3 10 


6 


500 


•099 


2 21 


12 


1,100 


•181 


2 9 


1 


2,100 


•285 


1 16 


2 


4,200 


•357 


1 


3 


6,300 


•379 


16 


Elastic 00 


GO 


•398 
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Spheroid with Eigidity of Iron (7*8 X 10^). 



Tjunar Semidiurnal llde. 




Modulus 




Reduction 


Acceleration 


of 


Yiscosity. 


of 


of 


relaxation. 




ocean tide. 


high water. 


Hrs. min. 






Hrs. min. 


Fluid 





•000 


3 6 


30 


140 


•420 


1 47 


1 


280 


•573 


1 7 


2 


560 


•647 


36 


3 


840 


•665 


25 


Elastic CO 


00 


•679 







Fortnightly Tide, 




Days, hrsp 






Days. hrs. 


Fluid 





'000 


3 10 


6 


1,700 


•294 


2 11 


12 


3,400 


•470 


1 18 


1 


6,700 


•602 


1 1 


2 


13,500 


•657 


13 


3 


20,200 


•669 


9 


Elastic 00 


00 


•679 






10. The influence of inertia, 



In establishing these results inertia has been neglected, and I will now show that 
this neglect is not such as to materially vitiate my results/'^' 

Suppose that the spheroid is Qor|.strained to execute such a vibration as it would do 
if it were a perfect fluid, and if the equilibrium theory of tides were true. Then 
the effective forces which are, according to D'Alembert's principle, the equivalent of 
inertia, are found by multiplying the acceleration of each particle by its mass. 

Inertia may then be safely neglected if the effective force on that particle which has 
the greatest amplitude of vibration is small compared with the tide-generating force 
on it. In the casQ of a viscous spheroid, the inertia will have considerably les^ effect 
than it would have in the supposed constrained oscillation. 

. Now suppose we have a tide-generating potential wr^^ cos (vt-^r)), then, according 
to the equilibrium theory of tides, the form of the surface is given by 



5 a 



2 



cr=^ B COS (vt+rj) ; 



^ In a future paper (read on December 19tli, 1878) I shall give an approximate solution of the problem, 
inclusive of the effects of inertia. 
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and this function gives the proposed constrained oscillation. It is clear that it is the 
particles at the surface whicb have the widest amplitude of oscillation. The effective 
force on a unit element at the surface is 

But the normal disturbing force at the surface is ^wa S cos {vt'^rj). Therefore inertia 

maybe neglected if ^ wv^ is small compared with 2wa, or if — t;^ is a small fraction. 

Ag 4:g 

The tide of the shortest period with which we have to deal is that in which v= 2a), so that 
we must consider the magnitude of the fraction 4 X -r— . If o) were the earth's true 

angular velocity, instead of its angular velocity relatively to the moon, then --: — would 

be the ellipticity of its surface if it were homogeneous. This ellipticity is, as is well 
known, -2^2-. Hence the fraction, which is the criterion of the negligeability of inertia, 

is about 5-8". 

If, then, it be considered that this way of looking at the subject certainly exag- 
gerates the influence of inertia, it is clear that the neglect of inertia is not such as to 
materially vitiate the results given above. 



II. 

A TIDAL YIELDING OF THE EARTH'S MASS, AND THE CANAL-THEORY 

OF TIDES. 

In the first part of this paper the equilibrium theory has been used for the determi- 
nation of the reduction of the height of tide, and the alteration of phase, due to bodily 
tides in the earth. 

Sir W. Thomson remarks, with reference to a supposed elastic yielding of the 
earth's body : " Imperfect as the comparisons between theory and observation as to 
the actual height of the tides has been hitherto, it is scarcely possible to believe that 
the height is in reality only two-fifths of what it would be if, as has been universally 
assumed in tidal theories, the earth were perfectly rigid. It seems, therefore, nearly 
certain, with no other evidence than is afforded by the tides, that the tidal effective 
rigidity of the earth must be greater than that of glass,"'^''* 

The equilibrium theory is quite fallacious in its explanation of the semidiurnal tide, 
but Sir W. Thomson is of opinion that it must give approximately correct results 
for tides of considerable period. It is therefore on the observed amount of the 
fortnightly tide that he places reliance in drawing the above conplusion. Under these 

* ^ Nat. PML', § 843, 
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circumstances, a dynamical investigation of the effects of a tidal yielding of the earth 
on a tide of short period, according to the canal theory, is likely to be interesting. 

The following investigation will be applicable either to the case of the earth's 
mass yielding through elasticity, plasticity, or viscosity ; it thus embraces Sir W. 
Thomson's hypothesis of elasticity, as well as mine of viscosity and elastico-viscosity. 

11. Semidiurnal tide in an equatorial canal on a yielding nucleus. 

I shall only consider the simple case of the moon moving uniformly in the equator, 
and raising tide waves in a narrow shallow equatorial canal of depth h. 

The potential of the tide-generating force, as far as concerns the present inquiry, is, 

- - sin^^cos 2((b'-'0)t), where t=-^ -r . This force will raise 

a bodily tide in the earth, whether it be elastic, plastic, or viscous. Suppose, then, 
that the greatest range of the bodily tide at the equator is 2E, and that it is retarded 

after the passage of the moon over the meridian by an angle-. Then the equation to 

the bounding surface of the solid earth, at the time ^, is r=:a+E sin^^ cos [2((^— co^)+€]; 
or with former notation or=Esin^^cos[2(</) — cot)'}-e]. 

The whole potential Y, at a point outside the nucleus, is the sum of the potential 
of the earth's attraction, and of the potential of the tide -generating force. Therefore 

Y=g"^+p(^J E sin^ 6 cos [2{<l.-cot) + e]+ ^^^' sin^ 6 cos 2(c^-o>0 

= c^^+ {F cos [2(.^-6)«)+e]+G sill [2(<^-ft)^)+e]}[^]%iii2 6 

3 T T . 

where F= -g^E^-- cos e, G= - sin e. 

Sir George Airy shows, in his article on ^' Tides and Waves" in the ^ Encyclopaedia 
Metropolitana/ that the motion of the tide-wave in a canal running round the earth is 
the same as though the canal were straight, and the earth at rest, whilst the disturb- 
ing body rotates round it. This simplification will be applicable here also. 

As before stated, the canal is supposed to be equatorial, and of depth h. 

After the canal has been developed, take the origin of rectangular coordinates in 
the undisturbed surface of the water, and measure x along the canal in the direction 
of the moon's motion, and y vertically downwards. 

We have now to transform the potential V, and the equation to the surface of 
the solid earth, so as to make them applicable to the supposed development. If v be 

the velocity of the tide-wave, then o)a=^v \ also the wave length is half the circum- 

2 
ference of the earth's equator, or 7ra ; and let m=-. Then we have the following 

transformations : — 
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Also in the small terms we may put r=a. Thus the potential becomes 

Y ^= const. -]-gy-{-F cos \j)i{x — 'z;^) + e]+Gsin[?7ite— i;^)+e]. 

Again, to find the equation to the bottom of the canal, we have to transform the 

equation 

7- = a+Esin^^cos[2(^— (i)^) + e]. 

If 2/' be the ordinate of the bottom of the canal, corresponding to the abscissa x, 
this equation becomes after development 

y'=zh — Ecos[m(x— i;i^) + €]. 

We now have to find the forced waves in a horizontal shallow canal, under the 
action of a potential V, whilst the bottom executes a simple harmonic motion. As 
the canal is shallow, the motion may be treated in the same way as Professor Stokes 
has treated the long waves in a shallow canal, of which the bottom is stationary. In 
this method it appears that the particles of water, which are at any time in a vertical 
column, remain so throughout the whole motion. 




Suppose, then, that x+^=x is the abscissa of a vertical line of particles PQ, 
which, when undisturbed, had an abscissa x. 

Let rj be the ordinate of the surface corresponding to the abscissa x. 

Let pq be a neighbouring line of particles, which when undisturbed were distant 
from PQ a small length k 

Conceive a slice of water cut off by planes through PQ, pg perpendicular to the 
length of the canal, of which the breadth is K Then the volume of this slice 

is 6 X PQ X Nn. 

Now PQ=/^— 'E cos [r}i{x—vt)-\-e]'—ri^ 

and Nn=^fl+-- 
Hence treating E and rj as small compared with h, the volume of the slice is 
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I CtdO 1% 1% J 



But this same slice^ in its undisturbed condition, had a volume hhh. Therefore the 
equation of continuity is 

dP 
^^zh-y—'E cos [m{x'—vt)-\-€]. 

Now the hydrodynamical equation of motion is approximately 

dY #i 



dx^' dx' dt^ 
The difference of the pressures on the two sides of the shce PQqp at any depth is 

d^ 

l^nX j^, ; and this only depends on the difference of the depressions of the wave- 
surface below the axis of x on the two sides of the slice, viz. at P and p. Thus 

df diq 

dx' ^ dx' 

Substituting then for j] from the equatioii of continuity, and observing that -^-^1 ^^ 

d^P 

very nearly the sanie as ^ |, we have as the equation of wave motion, 

fpt gy d^P 

gh~^^^-^mg E smlm{x'-vt) +e]= ^^ + ^ 

But 

dY 



dx' 
So that 



m F sin [w {x^ — t;^) + e] 4^ m G cos [m {x — ^ vt) + e]. 



j|==g/ij|+m{Gcos[m{i^'— t?i5) + e]— (F— E^) sin im{x-—vt)-^€]}. 

In obtaining the integral of this equation, we r^ay omit the terms which are 
independent of G, F, E, because they only indicate free waves, which may be 
supposed not to exist. 

The approximation will also be sufficiently close, \f ^ he written for x^ on, the right 
hand side. 

Assunae, then^ tl^at 

^"= A cos [r}%{x-^vt)-j-€]'^'B sin [m(a?— 1;^)+€]. 
By sub^titutioii in the equation of motion we find 

—m^{v^—'gh){A cos+B sin} =m{G cos— (1?— E^g) sin] . 
And as this must hold for all times and places, 

MDCCOLXXIX. E 
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G —l.aT sine 






B= 



F — 'Eg a(^^GO^e-~ f ^E) 



In the case of such seas as exist in the eai'th, the tide-wave travels faster than the 
free-wave, so that a^o}^ is greater than gh\ and the denominators of A. and B are 
positive. 

We have then — 

^r--.-——^ ^J (~ cose— -gE) sin— -sine COS y 

But the present object is to find the motion of the wave-snrfaoe relatively to the 
bottom of the canal, for this will give the tide relatively to the dry land. Now the 
height of the wave relatively to the bottom is 

PQz-y^— .E cos[m(a^— ^'^) + e]— '?; 



And 






clP 1 ( /t 

-f-^-r^ r-^ -cose' 

ax a-'oy^—gJi L \2 



~gE j cos + -z sin e sin > 



Hence reverting to the sphere, and putting a for a+hy we get as the equation to 
the relative spheroid of which the wave-surface in the equatorial canal forms part — 

Iv sm { T 2 

^.;=^_^-— ^l^cos 2(<^-cuO~g.gE cos [2(<^-a)^) + e] 

But according to the equilibrium theory, if V has the same form as above, viz. 

g^+^gl') Esin^(9cos[2(c^— <y?^) + e]+~ - sin^l9cos 2{(f>—o}t) 

and if r=a+'^ be the equation to the tidal spheroid, we have, as in Part I., 

sm \ T 3 

t^— — J -cos 2(^— cu?^)-f-r9'Ecos[2(<^— a)?^)-+-e] 

and the equation to the relative tidal spheroid is 

r=a+t^ — (T 

sin^ C T 2 

= a-| < -z cos 2{(j)—(ot) — -(^E cos [2{<f)—mt) + e] 

Now in either the case of the dynamical theory or of the eqnilibrium theory, if E be 
put equal to zero, we get the equations to the tidal spheroid on a rigid nucleus. A 
comparison, then, of the above equations shows at once that both the reduction of tide 
and the acceleration of phase are the same in one theory as in the other. But where the 
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one gives higli water, the other gives low water. The result is applicable to any kind 
of supposed yielding of the earth's nmass ; and in the special case of viscosity, the table 
of results for the fortnightly tide at the end of Part L is applicable. 

III. 
SUMMARY AND CONCLUSIONS. 

In § 1 an analogy is shown between problems about the state of strain of in- 
compressible elastic solids, and the flow of incompressible viscous fluids, when inertia is 
neglected ; so that the solutions of the one class of problems may be made applicable to 
the other. Sir W. Thomson's problem of the bodily tides of an elastic sphere is then 
adapted so as to give the bodily tides of a viscous spheroid. The adaptation is ren- 
dered somewhat complex by the necessity of introducing the effects of the mutual 
gravitation of the parts of the spheroid. 

The solution is only applicable where the disturbing potential is capable of expansion 
as a series of solid harmonics, and it appears that each harmonic term in the potential 
then acts as though all the others did not exist ; in consequence of this it is only 
necessary to consider a typical term in the potential. 

In § 3 an equation is found which gives the form of the free surface of the spheroid 
at any time, under the action of any disturbing potential, which satisfies the condition 
of expansibility. By putting the disturbing potential equal to zero, the law is found 
which governs the subsidence of inequalities on the surface of the spheroid, under the 
influence of mutual gravitation alone. If the form of the surface be expressed as a 
series of surface harmonics, it appears that any harmonic diminishes in geometrical 
progression as the time increases in arithmetical progression, and harmonics of higher 
orders subside much more slowly than those of lower orders. Common sense^ indeed, 
would tell us that wide-spread inequalities must subside much more quickly than 
wrinkles, but only analysis could give the law connecting the rapidity of the sub- 
sidence with the magnitude of the inequality.''^ 

* On this Lord Ratleigh remarks, tliat if we consider tlae problem in two dimensions, and imagine a 
number of parallel ridges, tbe distance between which is X, then inertia being neglected, the elements on 
which the time of subsidence depends are gw (force per unit mass due to weight), v the coefficient of 
viscosity, and X. Thus the time T must have the form 

T=(givyvyx\ 

The dimensions of giv, y, K are respectively ML~^T-^, ML-^T~\ L ; hence 

x-\-y=0 

V 

And aj= —1, y=l', ^= —1, so that T varies as • 

' *^ givX 

If we take the case on the sphere, then when i, the order of harmonics, is great, X compares with 7 ; 

i 

vi 

so that T varies as — -• 

ciwa 

E 2 



V 
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I hope at some future time to try whether it will not be possible to throw some 
light on the formation of parallel mountain chains and the direction of faults, by 
means of this equation. Probably the best way of doing this will be to transform the 
surface harmonics, which occur here, into Bessel's functions. 

In § 4 the rate is considered at which a spheroid would adjust itself to a new form of 
equilibrium, when its axis of rotation had separated from that of figure ; and the law is 
established which was assumed in a previous paper. '^''' 

In § 5 I pass to the case where the disturbing potential is a solid harmonic of the 
second degree, multiplied by a simple time harmonic. This is the case to be considered 
for the problem of a tidally distorted spheroid. A remarkably simple law is found 
connecting the viscosity, the height of tide, and the amount of lagging of tide ; it is 
shown that if v be the speed of the tide, and if tan € varies jointly as the coeffi- 
cient of viscosity aild v^ then the height of bodily tide is equal to that of the equi- 
librium tide of a perfectly fluid spheroid multiplied by cos e, and the tide lags by 

a time equal to 

It is then shown (| 6) that in the equilibrium theory the ocean tides on the yielding 
nucleus will be equal in height to the ocean tides on a rigid nucleus multiplied, by 

sin e, and that ther6 will be an acceleration of the time of high water equal to — — -. 

The tables in § 7 give the results of the application of the preceding theories to the 
lunar semidiurnal and fortnightly tides for various degrees of viscosity. A comparison 
of the numbers in the first columns with the viscosity of pitch at near the freezing 
temperature (viz., about 1*3 X 10^, as found by me), when it is hard, apparently solid 
and brittle, shows how enormously stiff* the earth must be to resist the tidally deform- 
ing influence of the moon. For unless the viscosity were very much larger than that 
of pitch, the viscous sphere would comport itself sensibly like a perfect fluid, and the 
ocean tides would be quite insignificant. It follows, therefore, that no very consider- 
able portion of the interior of the earth can even distantly approach the fluid state. 

This does not, however, seem to be conclusive against the existence of bodily tides in 
the earth of the kind here considered; for although (as remarked by Sir W. Thombof) 
a very great hydrostatic pressure probably has a tendency to impart rigidity to a 
substance, yet the very high temperature which must exist in the earth at a small 
depth would tend to induce a sort of viscosity — at least if we judge by the behaviour 
of materials at the earth's surface. 

In § 8 the theory of the tides of an imperfectly elastic spheroid is developed. The 
kind of imperfection of elasticity considered is where the forces requisite to maintain 
the body in any strained configuration diminish in geometrical progression as the time 
increases in arithmetical progression. There can be no doubt' that all bodies do 
possess an imperfection in their elasticity of this general nature, but the exact law 



* 



Phil. Trans., Vol. 167, Part L, sec. 5 of my paper. 
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here assumed lias not, as far as I am aware, any experimental justification; its 
adoption was rather due to mathematical necessities than to any other reason. 

It would, of course, have been much more interesting if it had been possible to 
represent more exactly the mechanical properties of solid matter. One of the most 
important of these is that form of resistance to relative displacement, to which the 
term ^^ plasticity " has been specially appropriated. This form of resistance is such that 
there is a change in the law of resistance to the relative motion of the parts, when the 
forces tending to cause flow have reached a certain definite intensity. This idea was 
founded, I believe, by MM. Thesca and St. Venant on a long course of experiments 
on the punchiog and squeezing of metals f"" and they speak of a solid being reduced to 
the state of fluidity by stresses of a given magnitude. This theory introduces a 
discontinuity, since it has to be determined what parts of the body are reduced to the 
state of fluidity and what are not. But apart from this difficulty, there is another 
one which is almost insuperable, in the fact that the differential equations of flow are 
non-linear. 

The hope of introducing this form of resistance must be abandoned, and the investi- 
gation must be confined to the inclusion of those two other continuous laws of resistance 
to relative displacement — elasticity and viscosity. 

As above stated, the law of elastico- viscosity assumed in this paper has not got an 
experimental foundation. Indeed, Kohlhausch's experiments on glass t show that 
the elasticity degrades rapidly at first, and that it tends to attain a final condition, 
from which it does not seem to vary for an almost indefinite time. But glass is one 
of the most perfectly elastic substances known, and, by the light of Thesca's experi- 
ments, it seems probable that experiments with lead would have brought out very 
different results. It seems, moreover, hardly reasonable to suppose that the materials 
of the earth possess much mechanical similarity with glass. Notwithstanding all 
these objections, I think, for my part, that the results of this investigation of the tides 
of an ideal elastico -viscous sphere are worthy of attention. 

There are two constants which determine the nature of this ideal solid : first, the 
coefficient of rigidity, at the instant immediately after the body has been placed in its 
strained configuration; and secondly, ''the modulus of the time of relaxation ot 
rigidity," which is the time in which the force requisite to retain the body in its 
strained configuration has fallen away to '368 of its initial value. 

In this section it is shown that the equations of flow of this incompressible elastico- 
viscous body have the same mathematical form as those for a purely viscous body ; so 
that the solutions already attained are easily adapted to the new hypothesis. 

The only case where the problem is completely worked out, is when the disturbing 

* " Sur recoulement des Corps Solides," Mem. des Savants Strangers, torn, xviii. and torn, xx., p. 75 
and p. 137. See also ' Comptes Eendus/ torn. 66, 68, and Liouville's Jonrn., 2^^*^ serie, xiii., p. 379, and 
xvi., p. 308, for papers on this subject. 

I POGGENDORFF Ann., voL 119, p. 337. 
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potential has the form appropriate to the tidal problem. The laws of reduction of 
bodily tide, of its lagging, of the reduction of ocean tide, and of its acceleration, are 
somewhat more complex than in the case of pure viscosity ; and the reader is referred 
to § 8 for the statement of those laws. It is also shown that by appropriate choice of 
the values of the two constants, the solutions may be either made to give the results 
of the problem for a purely viscoiis sphere, or for a purely elastic one. 

The tables give the results, for the semidiurnal and fortnightly tides, of this theory 
for spheroids which have the rigidity of glass or of iron — the two cases considered by 
Sir W. Thomson. As it is only possible to judge of the amoiint of bodily tide by the 
reduction of the ocean tide, I have not given the heights and retardations of the 
bodilv tide. 

It appears that if the time of relaxation of rigidity is about one quarter of the tidal 
period, then the reduction of ocean tide does not differ much from what it would be if 
the spheroid were perfectly ela,stic. The amount of tidal acceleration still, however, 
remains considerable. A like observation may be made with respect to the accelera- 
tion of tide in the case of pure viscosity approaching rigidity : and this leads me to 
think that one of the most promising ways of detecting such tides in the earth would 
be by the determination of the periods of maximum and minimum in a tide of long 
period, such as the fortnightly in a high latitude. 

In § 10 it is shown that the effects of inertia, which had been neglected in finding 
the laws of the tidal movements, cannot be such as to materially affect the accuracy of 

the results. 

[''" The hypothesis of a viscous or imperfectly elastic nature for the matter of the 
earth would be rendered extremely improbable, if the ellipticity of an equatorial 
section of the earth were not very small. An ellipsoidal figure with three unequal 
axes, even if theoretically one of equilibrium, could not continue to subsist very long, 
because it is a form of greater potential energy than the oblate spheroidal form, which 
is also a figure of equilibrium. 

Now, according to the results of geodesy, which until verj recently have been 
generally accepted as the most accurate — namely, those of Colonel A. R. Clarke t — 
there is a difference of 6,378 feet between the major and minor equatorial radii, and 
the meridian of the major axis is IS"" 34' E. of Greenwich. 

The heterogeneity of the earth would have to be very great to permit so large a 
deviation from the oblate spheroidal shape to be either permanent, or to subside with 
extreme slowness. But since this paper was read, Colonel Clarke has published a 
revision of his results, founded on new data; J and he now finds the difference between 
the equatorial radii to be only 1,524 feet, whilst the meridian of the greatest axis is 
8° 15' west. This exhibits a change of meridian of 24°, and a reduction of equatorial 

* Tlie part within brackets [ ] was added in JSTovember, 1878, in consequence of a conversation with 

Sir W. Thomson. 

t Quoted in Thomson and Tait, Nat. Phil., sec. 797. t Phil, Mag., August, 1878. 
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ellipticity to about one quarter of the formerly-received value. Moreover, the new 
value of the polar axis is about 1,000 feet larger than the old one. 

Colonel Clauke himself obviously regards the ellipsoidal form of the equator as 
doubtful. Thus there is at all events no proved result of geodesy opposed to the 
present hypothesis concerning the constitution of the earth. Sir W. Thomson 
remarks in a letter to me that ''we may look to further geodetic observations and 
revisals of such calculations as those of Colonel Clarke for verification or disproof of 
your viscous theory.''] 

In the first part of the paper the equilibrium theory is used in discussing the 
question of ocean tides ; in the second part I consider what would be the tides in a 
shallow equatorial canal running round the equator, if the nucleus yielded tidally at 
the same time. The reasons for undertaking this investigation are given at the 
beginning of that part. In § 11 it is shown that the height of tide relatively to 
the nucleus bears the same proportion to the height of tide on a rigid nucleus as in 
the equilibrium theory, and the alteration of phase is also the same ; but where the 
one theory gives high water the other gives low water. 

The chief practical result of this paper may be summed up by saying that it is 
strongly confirmatory of the view that the earth has a very great effective rigidity. 
But its chief value is that it forms a necessary first chapter to the investigation of the 
precession of imperfectly elastic spheroids, which will be considered in a future paper. '"^ 
I shall there, as I believe, be able to show, by an entirely different argument, that the 
bodily tides in the earth are probably exceedingly small at the present time; 



Appendix. 

November 7, 1878. 

On the observed height and phase of the fortnightly oceanic tide. 

In the following note I attempt to carry out the suggestion concerning the fort- 
nightly tide made in the preceding paper. 

The reports of the Tidal Committee of the British Association for 1872 and 1876 
contain the reductions of the tidal observations at a number of stations, into a series 
of harmonic tides, corresponding to the theoretical harmonic constituents of the tide- 
generating forces of the moon and sun. The tide with which we are here concerned 
is the fortnightly declinational tide. 

The heights of the tides at various times are all expressed in the form II cos {nt — e), 
where II is half the range of the tide in English feet, n the '' speed '' of the tide, and 
€ the retardation of phase, so that e-^7^ is the ^' lag'' of the tide. 

* Read before the Royal Society on December 19tli, 1878, 



32 MR. DARWIK OK THE BODILY TIDES OF VISCOUS AND SEMI-BLASTIO 

With the notation of the present paper n==20 for the fortnightly tide, and Ut is 
the '^ mean moon s ^' longitude from her node. 

The following are the results, giving the place of observation, its N. latitude, and 
the years of observations. With respect to Brest and Toulon, R is reduced to feet 
from centimetres, so as to be naade comparable with the other results : — ■ 



Ramsgale, 
about 5 r 21^. 


Liverpool, 
53° 40^ 


Hartlepool, 

54° 41^ 


1864. 


1857-58. 


1858-59. 


1859-60. 


1866-67. 


1858-59. 


1859-60. 


1860-61. 


R ^0331 

€ 268°-29 


•093 

170°^7 


•037 
148°'8 


'024 
72°-9 


•036 
340° • 6 


•052 
190"^34 


•053 

222° -34 


•073 
158°-62 



Brest, 
48^ 23'. 


Toulon, 

43° 7'. 


Kurrach^e, 
24° 53'. 


Cat Island, 

Gulf of Mexico, 

80° 23'. 


1875. 


1853. 


1868-69. 


1869-70. 


1870 71. 




R ^099 

6 80°^65 


•051 
.139°-50 


•038 
335°-40 


•064 
333" -91 


•035 

283° ^22 


•043 
136^-69 



Tn their present form the observations do not appear to present any semblance 
of law, but when they are rearranged we shall be able to form some idea as to whether 
they are really quite valueless ox not for the point under consideration. 

The theoretical expression for the fortnightly tide of 'm ocean covering the whole 
earth, according to the equilibrium theory, is 



3 T 



Wnere 



,^-a sin^ i(i-—co,s^ 6) cos 2Xlf, 
10 2 \6 / 

~^, Ct=-^, a= earth's radius, i the averas^e obliquity of the earth's axis to 
2c'^'^ ba ^ ^ ^ 



the normal to the plane of the lunar orbit during the fortnight in question^ 6 the 
colatitude of the place of observation, 

If we take ^=23'^ 28' the obliquity of the ecliptic, a— 20 "9 million feet, we find 



3 T 



10 S 



asin^r='207 foot, 



So that the fortnightly tide should be expressible by 

•207 (i— sin^ (lat)) cos 2m, 
In Thomson's corrected equilibrium theory the second factor should be 

X "t" ^^ • 0/1 1 \ 

— --— — sm^ (lat,) 
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where ® is a certain definite integral, depending on tlie distribution of land and 
water^ but wHch has not yet been evaluated. 

The latitude of evanescent fortnightly tide is 36° 15' if (S is zero ; and if we bear 
in mind that ® may be negative, it is clear that the observations at Cat Island 
(lat. 30° 23') are made too near the critical latitude to be trustworthy for determining 
the true fortnightly tide. It is also hardly possible to believe that the observations at 
Tou.lon should show a true tide of this denomination, because the Mediterranean must 
be regarded as a virtually closed sea. The observations at Cat Island, and at Toulon, 
will therefore be set aside. 

The first process to be applied to the above observations is obviously to divide each 
value of R by ^ -^ sin^ (lat.) ; the following are the factors for reducing the values 
ofE:— 



[i^sin^Clat.)] 



-1 



Hamsgate. 


3-62 


Liverpool. 


3-17 


Hartlepool. 


3-01 


Brest. 


3-07 


TCurracliee. 


6-40 



These factors will be applied to the values of R in the table first given. 

The next point to consider is the phase of the tide. The formula we have given 
shows that the fortnightly tide consists in an alternate deformation of the ocean level 
into an oblate and prolate spheroid of revolution, when the tide is deemed to be 
superposed on a true sphere, instead of on an oblate nucleus. 

When t is zero the spheroid is oblate, and this may be called high-tide; when t-=:^~- 

it may be called low-tide. It follows, therefore, that N. of lat. 36° 15' high-i(^(fe is 
\o'^ -water ^ and vice-versd; but S. of this latitude the tide and water agree. But the 
formulae in the tidal reductions always refer to Ingh-iuater, hence to find the retarda- 
tion of the tide we must subtract 180° from all the e's for places N. of 36° 15' — that is 
to say (Cat Island being rejected) for all except Kurrachee. 

For Kurrachee, we may observe that any retardation e may be regarded as a retarda- 
tion e — Stt, which, if negative, is an acceleration of tide. If 27r — e be less than 180°, 
this appears to be the more correct light in which to look at it. 

Now if we reduce all the observations in the way indicated, so that the fortnightly 
tide is given by Il'(-|— cos^ 6) cos (20^—07), we find the following results : — 
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Ramsgate. 


Liverpool. 






Hartlepool. 




R' -120 
V +88°'29 


•295 
9°-3 


•117 
3r'2 


•076 
107^-1 


•114 

+ 160"-6 


•157 

4-10°-34 


•160 

+ 42°- 34 


•220 
2r-38 
















Brest. 


Kurrachee, 






R' ^304 
7) 99°-35 


•243 
24° •GO 


•410 

1 26°^09 


•224 
76°^78 





We will consider R' first. 

From these twelve values we find R'='203^ with a probable error ±'068. 

The value of R'is almost exactly that indicated by theory (viz., '207), but the very 
large probable error renders the result so uncertain, that it can only be asserted that 
the results do not disprove a diminution of fortnightly tide. 

With regard to phase, it will be observed that there are eight cases of accelerated 
tide to four of retarded. Two of the retarded tides refer to Hartlepool, and concerning 
this station Sir W. Thomson" says in the report : ^' There is scarcely suflficient agree- 
ment between the results deduced from the long-period tides to be satisfactory, 
although the quantities of some are within reasonable limits.'^ 

It maybe remarked, in passing, that Cat Island gives a retarded tide, and Toulon an 
accelerated one. 

If we treat these alterations of phase in the same way as R^ was treated, we find a 
mean acceleration of phase of 7^*85, but with a p.e. several times larger than the 
result itself. But, in fact, with so few and such irregular observations the method 
of least squares is useless. 

The cases of retarded phase certainly show considerably more irregularity than 
those of accelerated phase. If we take the mean of the cases with accelerated phase, 
we shall find an acceleration of 48^, which corresponds in time to an acceleration 
of 1 day 20 hours. 

Now three out of four years of observations show an accelerated tide at Liverpool ; 
all three years show an acceleration at Kurrachee ; the Hartlepool observations are 
not of very much value ; while the single year of acceleration at Brest may be set off 
against the single year of retardation at Bamsgate. If then we ask ourselves whether 
acceleration or retardation is the more probable, I think it must be answered in favour 
of acceleration ; and if so there seem to be some indications of a viscous yielding of 
the earth's mass. It must be admitted, however, that the evidence is exceedingly 
uncertain. 

It does not seem to be noticed in the tidal reports, that amongst the '^ Helmholtz 
compound shallow- water tides " there will be found several which have the same 
period, or very nearly the same period, as the true fortnightly declinational tide. If 
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we write (as in the report) y, a for earth's rotaton and moon s mean motion, then we 
shall find the following speeds will combine so as to give shallow-water tides indis- 
tinguishable from the true fortnightly tide, namely, 2(y— cr) and 2y, also y— 2cr and y ; 
and besides there are four combinations of the elliptic tides which give the samq 
period of compound tide, if we neglect the motion of the moon s perigee. It therefore 
seems quite possible that in certain stations the true fortnightly tide may be masked, 
or have its phase largely affected by these compound tides, and this is, perhaps, the 
explanation of the great irregularity in the phases. 

A series of observations in some oceanic island near the Equator, or better still fa,r 
north or south, would be of immense value to decide this point. 



F 2 



